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Metal atoms adsorbed on few layer graphenes condense to form nanometer-size droplets whose 
growth is size limited by a competition between the surface tension and repulsive electrostatic in- 
teractions from charge transfer between the metal droplet and the graphene. For situations where 
the work function mismatch is large and the droplet surface tension is small, a growing droplet 
can be unstable to a family of shape instabilities. We observe this phenomenon for Yb deposited 
and annealed on few layer graphenes and develop a theoretical model to describe it by studying 
the renormalization of the line tension of a two dimensional droplet by repulsive interparticle in- 
teractions. Our model describes the onset of shape instabilities for nanoparticles where the growth 
is size-limited by a generic repulsive potential and provides a good account of the experimentally 
observed structures for Yb on graphene. 

PACS numbers: 68.65.Pq, 68.70.+W, 68.37.Lp, 68.43.Hn 
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I. INTRODUCTION 

Graphite is widely used as a substrate for the synthesis 
of free standing metal nanoparticles due to its chemical 
inertness, low diffusion barriers for adsorbed species and 
its compatibility with various electron microscopies. In 
earlier work [l| we observed that when Au atoms are 
deposited on few layer graphenes (FLG) containing m- 
layers where 1 < m < 20, the mode of adatom con- 
densation and nanoparticle growth differs sharply from 
that observed on thick graphite. On few layer graphenes 
Au condenses to form size-limited isotropic droplets as 
shown in Fig. [T] where the diameter is controlled by the 
layer count m of the graphene substrate, increasing ap- 
proximately as m}'^ . This phenomenon is well described 
by a model in which electrostatic dipole-dipole repulsion 
within a condensed island prevents the continued growth 
of large droplets. The thickness dependence arises from 
a short range cutoff of this repulsive potential which oc- 
curs on the scale of the dipole. For few layer graphenes, 
which screen poorly when the charge exchange is small, 
this scale is effectively the thickness of the graphene film. 



In this paper we examine the effects of general re- 
pulsive interactions on shape selection for a condensed 
two dimensional system. This work is motivated by ex- 
periments similar to [l|] where Yb is adsorbed on few 
layer graphene films. As shown in Fig. [5] these are 
observed to condense in filamentary labyrinthine struc- 



tures rather than in isotropic size-limited droplets, sug- 
gesting new physics in the surface energetics. Indeed, 
Yb is distinguished by a significantly larger workfunction 
mismatch to the graphene and a lower intrinsic surface 
tension. These two features produce a mode-dependent 
renormalization of the line tension of a growing isotropic 
droplet. A shape instability occurs when the line ten- 
sion for a particular surface mode of the droplet, con- 
strained to have a constant total area, goes to zero. We 
develop a new formulation describing this shape insta- 
bility driven by a general repulsive potential. For ap- 
plication to Yb/graphene, where this repulsive potential 
repulsion is provided by electrostatic dipolar interactions 
and can be calculated from the work function mismatch, 
we observe a sign reversal of the mode dependent surface 
tension for an isotropic droplet as a function of its radius. 
This model correctly accounts to the observed widths of 
the labyrinthine structures and the occasional appear- 
ance of threefold and fourfold vertices in the condensed 
filamentary patterns. 

Shape instabilities produced by electrostatic and mag- 
netostatic dipolar interactions often arise in condensed 
matter problems. For example, in two dimensions they 
are associated with the shapes of amphiphile domains at 
the air water interface (electrostatic) [2| and of ferrofluid 
droplets (magnetostatic) j3|. Langer et al. [4| examined 
this problem theoretically, exploiting the dipolar form of 
the repulsive interaction. They derived an explicit for- 
mula for the line energy of a two dimensional droplet in 




FIG. 1: SEM image of Au islands formed after annealing 
on four layer graphene. The Au atoms condense to form size 
limited nearly isotropic droplets whose radius is limited by the 
repulsion between perpendicular dipoles at the Au graphene 
interface. The average droplet radius is determined by the 
microscopic width of the dipole layer, effectively the width of 
the few layer graphene substrate. 



terms of a double line integral over its perimeter, a result 
that in principle can be used as input to an algorithm to 
compute the shape of a growing droplet. This approach 
has been refined by Iwamoto and coUaborators [5|, ISj who 
replace the double line integral by an expression for the 
hne energy parameterized by the amplitudes of its modes 
of deformation, and they applied this method to study 
effects of dipolar interactions both for the perpendicular 
geometry and for dipoles tilted with a nonzero component 
parallel to the tangent plane of the droplet. Our work is 
similar in spirit to that of Ref. [5[, though we present 
it here in a new form that is applicable to a general 
repulsive potential. For the electrostatic dipolar repul- 
sive potential appropriate to the Yb/graphene problem 
the model provides a good account of the experimental 
observations. Furthermore this formulation emphasizes 
that shape instabilities of this type are a generic prop- 
erty of condensed phases of species with repulsive tails in 
their interaction potentials. Our method can be used to 
access this physics for a droplet with a general two point 
interaction potential in its interior. 

Section II of this paper provides more information 
about our sample preparation, characterization and 
imaging of Yb and Au nanoparticles formed on few-layer 
graphenes. Section III briefly reviews the isotropic model 




FIG. 2: SEM image of Yb nanoparticles formed after an- 
nealing on graphene. The metal atoms condense to form 
anisotropic fibrillar structures. The top panel is for an av- 
erage coverage of 0.12 nm and the bottom is for an average 
coverage of 0.5 nm. 



introduced in [l| appropriate to size-limited circular 
droplets for Au on graphene. Section IV presents some 
useful formulas for weakly perturbed circular droplets. 
Section V constructs a model for the droplet energet- 
ics, including the interaction renormalized line tension. 
Section VI applied the model to Yb and Au droplets on 
few layer graphene and provides comparison of the model 



with the experimental data. A brief discussion of the re- 
sults in given in Section VII. 



II. EXPERIMENTAL 



We prepared graphene flakes by mechanical exfolia- 
tion of kish graphite onto 300 nm silicon oxide on silicon 
wafers. We then cleaned the sample by annealing at 400 
C in a reducing atmosphere, 1:1 H2 and Ar. Individ- 
ual flakes were identified by color contrast in an optical 
microscope. Flake thickness was determined by Veeco 
Dimension 3100 AFM in tapping mode. 

We then thermally evaporated thin layers of metal onto 
the surfaces. The deposition thickness was determined 
by a crystal thickness monitor. We deposited 0.22 ± 0.1 
nm and 0.5 ± 0.1 nm of Yb on different chips. For Au 
samples, the thickness was 0.3±0.1 nm. Each sample was 
annealed to equilibrium at 600 C (Yb) or 400 C (Au) for 
three hours in a reducing atmosphere, 1:1 H2 and Ar. 

We then imaged the samples in an FEI DB 235 high- 
resolution SEM / FIB in HRSEM mode. Even short 
exposures to the beam have been found to contaminate 
the imaged region and reduce quality, so we were espe- 
cially careful that each image was taken in a previously 
unimaged area. 

Figures 1 and 2 show Au and Yb nanoparticles on 
single-layer and few layer graphenes. The Au particles, 
as noted in [1], are nearly isotropic. The average radius of 
these nanoparticles increases with the thickness of the few 
layer graphene substrate. By contrast Yb forms branch- 
ing strands. The sample with less Yb has strands that 
are more regular in width. 

Figure 3 shows histograms of the distribution of Yb 
strand widths extracted from the images of Figure [2j The 
data at higher coverage were obtained over a wider image 
than shown in the lower panel of Figure [2] We excluded 
only places where the Yb was clearly branching. The 
mean of the distribution for 0.22 nm average coverage is 
4.01 ± 0.08 nm, and at a 0.5 nm coverage it increases to 
is 4.95 ±0.11 nm. The peaks of the two distributions 
overlap; the increase arises from the irregularities which 
are more marked on the sample with more Yb. We also 
measured the orientations of those Yb strands that have a 
clear direction. There is no clear deviation from isotropy 
on the whole. 
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FIG. 3: Width distribution of the Yb fibrils displayed in the 
top and bottom panels of Figure [2] for the average coverages 
shown. 



(Figure S]) can be expressed 

C/ = 7 {nR^ + 27r i?/i) -|- CTrR^h + Fi?^ 

where 7 is the surface tension and e (< 0) is bulk cohe- 
sive energy density. The last term results from the dipole- 
dipole interactions and the scaling with _R^ is the essential 
feature that prevents the growth of large droplets. Note 
that because of the long range tail of the electrostatic 
interaction, the last term can be regarded as a Volume^ 
contribution to the droplet energy. Evaluation of the 
coefficient F in this expression requires that one regular- 
ize the short range divergence of the l/r^ dipole-dipole 
potential. Physically this is regularized by the finite spa- 
tial extent of the microscopic dipole (which we label d). 
For graphene doped near its charge neutrality point the 
screening is weak and the distance d may be identified 
with the thickness of the graphene film [7[. This leads 
to the central result that the droplet radius is propor- 
tional to d^/"^ which is the scaling rule identified in the 
experiments studying Au particles on graphene [l| . 

The electrostatic contribution in this expression is ob- 
tained in [ll by isolating the small momentum {q -^ 0) 
limit of the dipole-dipole interaction. By contrast, the 
model for the shape instability developed below will re- 
quire us to extend this into the q ^ regime, and in 
fact the interactions that drive the shape instability are 
obtained by integrating the interaction over all momenta. 



IV. MENSURATION FORMULAE 



III. REVIEW OF THE ISOTROPIC MODEL 

Reference [l| considers the contribution to the energy 
of an isotropic droplet, modelled as a short circular cylin- 
der. The energy of a droplet of radius R and height h 



We are interested in perturbations of the droplet 
around a reference circular shape. Referring to Fig. |4l we 
write the droplet radius as a function of the polar angle 
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Any shape change of a droplet that preserves the parti- 
cle number keeps the norm A constant (assuming fixed 
h). In particular, if we consider an area-preserving de- 
formation characterized by the amplitudes r^ then the 
isotropic term has to adjust in the manner r^ = ro[r„i] = 
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FIG. 4: Diagram illustrating the geometry of the deformed 
droplet with radius r{(j>), height h and dipole layer height d. 
The lower panel illustrates the three modes of deformation 
with m = 2,3,4. 



Since r{(f>) is real r_m — r^. We are considering insta- 
bilities around the circular shape and so we consider the 
situation r^ <C tq. 

The differential length along the tangent line of this 
curve is 



de-" 



dr^ + {rd(ff 



which gives a formula for the length of the edge of the 
droplet 
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V. DROPLET ENERGY 



V.l Line Energy 



The effective line tension T = jh giving for the line 
energy 

C/iine = 2TT"fh r(, + — ^ m^jr^p J 

where we are truncating the expansion at quadratic order 
in the r™. The constant area constraint can be enforced 
by writing 



C/linc = 2TTjh '"O + ^ ^ ("1^ - l)kn 



Note that to quadratic order in the deformations, the 
m — 1 term does not change the line energy since 
it describes a rigid translation of the reference circular 
droplet. 



V.2 Interaction Energy 



Note that because of the </> derivative the higher m terms 
are more effective at increasing the arclength and thus 
tend to be suppressed by a positive surface tension. 

The cross sectional area is given by a line integral over 
the perimeter 
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Since r[(f) — r{(j})er one has 



which gives 



,-, dr ^ 

dr — -r-Br 
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r(0)2 d0 



and integrating over angles gives 

A = Tirl + TT 2J kri 



Two-point interactions inside the condensed droplet 
can modify the effective line tension at its boundary. In 
this section we develop a formalism for calculating the 
mode-dependent renormalized line tension working in a 
constant-area ensemble. One can also derive these results 
by working in the grand canonical ensemble, allowing for 
fluctuation in the total droplet area. 

We describe the interactions within the droplet by a 
two dimensional potential which depends on the lateral 
separation of two particles V{ri — r2)- In this prob- 
lem we are particularly interested in the interactions be- 
tween electric dipoles mutually oriented along the inter- 
face normal with dipole density r. This dipole layer pro- 
duces a potential step that equilibrates the work func- 
tion mismatch between the metal and graphene A$ with 
T — A<I>/47re. The electrostatic energy of an island can 
be expressed as a double integral over the area of the 
droplet 



Ud = 



A$2 
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rfVdV/(f)/(r')5(|f-r'|) 



where the dimensionless distribution function /(r) — 1 
inside the droplet and zero outside. It is useful to write 
this interaction energy as a momentum integral 



where 



and 
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At distances large compared to the size of the dipole d, 
the interaction kernel takes the form 

lim g(r) -^ -^ 

which must be regularized in the near field on the scale of 
the dipole. A convenient form for the regularized kernel 
is 



ffW 



1 



(r2 + d2)3/2 



which gives 



v{q) = ^ cxp(-qd) = -w(^) 



where S, — qd defines the dimensionless momentum. Note 
that the interaction strength scales inversely with the 
width of the dipole layer d. 

The function S{q) can be parameterized in terms of 
the coefficients r„i. Consider the Fourier integral over 
the area of the droplet 
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Here the isotropic part gives 



e-"'-'' r dr d<j) 



Saiq) = 27r / Jo{qr)rdr 
Jo 
27rro 
= Jiiqro) 

while the anisotropic piece can be decomposed by noting 
that rm ■C ^o and thus 

5„ = 27rz"e™'^'J„(qro)ror„ 



and therefore we have that 
S{q) = ^^Ji(gro) + 27r^z'"e™'^'J™(<7ro)ror„ 

m 

Using this in our expression for the interaction energy, 
we find that it can be partitioned into an isotropic part 
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and a shape dependent piece which is given by a sum 
over m 
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We observe that the coefficient of the 



term is always 



positive and thus C/ghape is positive definite. A shape in- 
stability arises from a competition between this energy 
and compensating changes to the isotropic interaction 
energy which necessarily occur if the droplet area (equiv- 
alently the number of condensed particles) is held fixed. 



V.3 Area Preserving Deformations 

We wish to isolate from the interaction energy the part 
that can be associated with a simple renormalization of 
the surface tension (it depends only on the droplet area) 
and a residual part which contributes to the effective line 
tension (depends on the droplet shape). To do this we 
write 

U=U± JintA =tj + 7int^ 

where A = 7r(rQ + X^mkmP)- ^^ appropriate choice 
of 7int removes from the energy difference U its area- 
dependent part and isolates the pure shape dependent 
energy. U can be expressed in terms of 7int and the 
coefficients [gm] of the quadratic terms in the interaction 
energy 

U = (go - Trjint)rl + ^ (.g™ - 7r7i„t)|rmp 



Setting dU/dA = eliminates the differential area de- 
pendence from U and solving for 7int then gives the con- 
tribution of the interaction term to the area-dependent 
energy. This gives 



dU 
dA 
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7int = 




The interaction contribution to the surface tension is 



7int 



The residual energy U = J^m (dm - 7r7i„t) \r^\'^ contains 
the shape-dependent energy under the constraint of con- 
stant total droplet area. 

In the evaluation of 7int care must be taken to retain 
all terms that contribute to the energy U at order |rmp, 
noting that the expansion coefficients [gm] are themselves 
function of rp. The 77i-dependent terms in 7int^ can then 
be organized in ascending powers of the aspect ratio ro/d 
in the manner 

ro\ JH^ro/d) 



e 



(7)' "^i(^^o/^)( 



M^rp/d) - J2{^rold) 
2^ 



Vd) ^^'"(^'^0/^''*- 2 ' 

When ro/d is large, this expression is controlled by its 
terms containing the highest powers of r^/d. We ob- 
serve that the first term at ©(rp/d)^ exactly cancels the 
"bare" contribution of the ?7i-th mode to the interaction 
energy. The remaining term at O(ro/(i)^ ultimately dom- 
inates the momentum integral, since the last term is a 
rapidly oscillating function of f as illustrated in Fig. [S] 
The momentum-integrated interaction is seen to be neg- 
ative indicating that for area-preserving deformations re- 
pulsive interparticle interactions destabilize the isotropic 
shape at sufficiently large radius. To leading order in 
the small quantity d/ro we find that the interaction term 
simplifies to 



U, 



shape 
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V.4 Renormalized Line Tension 



Combining our results for the bare and interaction 
renormalized line energies, for a droplet constrained to 




FIG. 5: Momentum dependence of the integrand for evalua- 
tion of shape dependent interaction energy for m = 2 plotted 
as a function of the dimensionless momentum ^ = qd, broken 
into separate contributions /p sorted by their dependence on 
{ro/dy (top panel) and the total (bottom panel). The thin 
(red) line in the bottom panel plots an expression discussed in 
the text that provides a good approximation to the integrated 
coupling strength. The data are presented for ro/d = 5. 



have constant area we have 



Celine = ^ TTm^jh- 



m^O 



ro 
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where z = ro/d and we define the dimensionless coupling 
constant 



A: 



A$2 



and the weight function for the TTi-th mode 

31 + M^Z)J2{^Z) 



Fra{C,z) = -z' 



A^Jm{^z){Jm-l{£,z) " Jm+l{^z)) 



Because of the rapid oscillation of the Jm's for large rg/d 
the integrated coupling strength is well approximated by 
retaining only the constant term in the z^ coefficient; in 
this approximation the momentum-integrated interaction 
is independent of m and we have 



A. 



1 
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i + 4f^^' 



The lower panel of Figure [5] compares the complete inte- 
grand of the momentum integral to an approximate form 
that retains only the constant in the z'^ term. Vanishing 
Am signifies the onset of a shape instability in the m-th 
deformation mode. Thus the critical radius for destabi- 
lizing the TO-th mode is 



rc{m) f2{rn? — 1) 



1/3 




VI. COMPARISON WITH EXPERIMENT 

For Yb on graphene we have A<i> — 2.2 eV [8|, 7 ~ 
320 erg/cm^ [^ and h w 10"'^ cm which gives A sa .067. 
Figureiniplots the renormalization coefficients Am for this 
coupling strength as a function of rg/d for m — 2,3 
and 4, computed using the full expression for 7inf The 
plot shows that Yb islands grows from an isotropic seed 
to a radius ro/d ^ 4.5 where the quadrupolar m ~ 2 
mode becomes unstable. (The approximate expression 
gives rc{m = 2)/d k, 4.47 and is indeed very accurate 
in this regime.) Above this critical point one expects 
an exponential growth of the droplet along a single axes, 
which one can associate with the filamentary structures 
observed experimentally. The width of such a filament 
is twice this critical radius; for d sa 0.5 nm we find that 
filament width is ~ 5 nm in very good correspondence 
with the average width 4.01 ± .08 nm computed from the 
distribution shown in Fig. |3]at 0.12 nm average coverage. 
As the coverage increases the filaments fold and coarsen 
due mode competition between the fibrils. The average 
width in the higher coverage state is thus slightly larger 
(4.95±0.11 nm). One can contrast this with the situation 
for Au adsorption ongraphene [l| (where A4> w 0.5 eV, 
7 « 1130 erg/cm^ [10|, and A ^^ xlO~^) which gives 
Tcim — '2,)/d K, 20. This requires growth of a circu- 
lar droplet to a diameter exceeding 20 nm which is larger 
than both the intrinsic size limit imposed by the isotropic 
term in the dipolar energy and the largest droplet sizes 
observed experimentally for Au on few layer graphenes 

a. 

In general the stability limits for the isotropic and dis- 
torted drops show different scaling with the width of the 
dipole layer d. A size-limited circular droplet has a ra- 
dius that increases oc d^/^ while it is stable against shape 
fluctuations below a critical radius oc d. Thus in the limit 



FIG. 6: Line tension renormalization coefficients Am plotted 
as a function of ro/d for m — 2,3,4: and for A « .067, appro- 
priate to Yb/graphene. 



of weak coupling (large d) a size-limited circular droplet 
is stable, while for strong electrostatic coupling with d ~ 
interlayer spacing the droplet can undergo an interaction- 
driven shape instability. Au and Yb provide examples, 
respectively, of this weak coupling and strong coupling 
behavior. 



VII. DISCUSSION 

Three features of Fig. [6] are noteworthy. First, the in- 
stability is driven by the terms in our expansion for the 
shape energy that are quadratic in the deformation am- 
plitudes and have the strongest (highest power) depen- 
dence on ro/d. Thus the instability requires only growth 
of a droplet past a critical radius for the shape change 
to occur. At larger r^/d there is no mechanism that can 
restore a stable isotropic solution. Thus sufficiently large 
droplets are absolutely unstable to this type of shape in- 
stability. Second, the renormalization coefficient shows 
that the crossover from the weakly renormalized regime 
(small To) regime to the unstable regime (large ro) oc- 
curs over a narrow size range. This also reflects the very 
strong z^ dependence of the dominant term in the inte- 
grated coupling strength. Thus one expects the fibrils to 
show a sharply peaked width distribution, as is demon- 
strated experimentally in Fig. [3l Third, it is striking 
that the integrated coupling strength is nearly the same 
for all the modes of deformation of the circular droplet. 
Ultimately the shape instability is suppressed for large m 



modes because of the m? scaling of the bare hne tension 
rather than through the residual m dependence in the in- 
teraction contribution. This leads to our simple scaling 
rule for the m-dependent critical radii. 

Finally, we note that previous continuum formulations 
of this problem inevitably require a finite droplet height 
to regularize the short distance singularity in the dipole- 
dipole potential when they are treated as point dipoles 
[J, |5|. In our treatment this is regularized more nat- 
urally by representing the interaction potential using a 
nonsingular near field form controlled by the finite size 
of the relevant microscopic dipoles. More generally our 
expression for the interaction-renormalized line energy 
now can be applied to any repulsive two point potential 
in the droplet interior. Indeed the scaling form for F^ 
describes the renormalization of the line tension by any 
nonsingular bulk interaction. 
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